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摘 要 :微细 铣 涡 中 ,由 于 受 加 工 余 量 限 制 , 铣 刀 每 元 进 给 量 非常 小 甚至 小 于 材料 内 部 晶 粒 尺寸 。 
当 切 削 在 晶 料 间 进 行 时 会 造成 切削 不 连续 ,从 而 使 铣削 过 程 颤 振 和 急剧 增 大 ,引起 加 工 表面 质量 下 
降 ,严重 影响 微型 零件 的 疲劳 寿命 和 使 用 性 能 。 本 研究 基于 Voronoi 图 建立 多 晶体 几何 模型 ,通过 
仿真 与 实验 验证 的 方法 对 刀具 刃 口 圆 陶 半径 、 切 前 深度 和 蝇 粒 尺寸 对 已 加 工 表面 残余 应 力 的 影响 
规律 进行 详细 研究 ,建立 多 元 线性 回归 模型 并 进行 显著 性 检验 。 对 经 热处理 后 具备 不 同 唱 粒 尺 寸 
的 高 温 合 金 GH4169 进行 微 铣削 实验 ,通过 X 射线 衍射 法 测量 工件 进 给 方向 已 加 工 表面 的 残余 应 
力 值 ,将 实验 值 与 预测 值 进行 对 比 ,从 而 验证 残余 应 力 预 测 模型 的 准确 性 。 结 果 表 明 , 其 误差 分 布 
在 2.18% ~7.35% ,不 超过 10% ,说 明 残 余 应 力 的 预测 模型 准确 性 较 高 。 
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Research on residual stress in micro-milling of GH4169 
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Abstract:In micro-milling, due to the limitation of machining allowance ,the feed per tooth of the milling 
cutter is very small or even smaller than the internal grain size of the material. When the cutting is carried 
out between the grains ,it will cause discontinuity in cutting ,which will sharply increase the chatter during 
the milling process ,cause the quality of the processed surface to decrease ,and seriously affect the fatigue 
life and performance of the micro-parts. Based on the Voronoi diagram ,the paper establishes a polycrystal- 
line geometric model. Through the method of simulation and experimental verification ,the influence law of 
the cutting edge arc radius ,cutting depth and grain size on the residual stress of the machined surface is 


studied in detail ,and a multiple linear regression model is established. Perform a significance test. The 
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micro-milling experiment was carried out on the superalloy GH4169 with different grain sizes after heat 


treatment. The residual stress value of the machined surface in the feed direction of the workpiece was meas- 


ured by X-ray diffraction method ,and the experimental value was compared with the predicted value to veri- 


fy the residual stress prediction model Accuracy. The results show that the error distribution is 2. 18%- 


7. 35% ,not more than 10% ,which can indicate that the prediction model of residual stress is more accurate. 


Key words :grain size;micro-milling;residual stress; prediction model 


在 微细 铣削 中 切削 参数 减 小 ,与 材料 的 晶 粒 尺 
才 相 当 , 实 际 上 的 切削 相当 于 在 不 同 晶 粒 中 进行 , 造 
成 切削 过 程 的 不 连续 ,工件 应 该 被 视 为 不 均匀 的 、 各 
向 异性 的 。 微 型 零件 由 于 尺寸 小 、 刚 度 低 , 加 工 过 程 
中 容易 产生 变形 ,从 而 导致 其 可 靠 性 及 使 用 性 能 降 
低 0 2 。 微 型 零件 的 使 用 性 能 与 工件 的 表面 质量 息 
息 相 关 , 在 表面 质量 的 各 个 指标 中 ,已 加 工 表面 残余 
应 力 对 零件 的 静 力 强度 , 抗 疲劳 强度 ,使 用 寿命 等 有 
着 重大 影响 ?] 。 在 对 微细 铣削 过 程 研究 中 ,考虑 工件 
材 狂 的 晶 粒度 ,对 切削 力 的 控制 和 切削 完成 后 零件 的 
表面 质量 控制 都 有 着 重要 的 指导 意义 ,建立 多 晶体 材 
料 杖 切 前 模型 可 以 更 准确 地 预测 最 小 切 悄 厚度、 切削 
太 我 余 应 力 等 ,对 实际 加 工 过 程 具有 指导 作用 。 
CD 传统 切 前 中 ,切削 参数 相对 较 大 ,可 以 忽略 晶体 
续 梅 中 的 点 缺陷 .位 错 缺 陷 和 裂纹 等 微观 缺陷 ,工件 
下 以 被 视 为 均匀 的 、 各 向 同性 的 ,多 晶体 材料 由 于 塑 
性 变形 导致 颤 振 急剧 增 大 而 影响 加 工 过 程 中 切削 力 
和 现 削 温度 ,从 而 对 残余 应 力 的 影响 不 可 以 忽略 。 
通过 研究 刀具 刃 口 圆 弧 半 径 与 材料 晶 粒 尺寸 对 切削 
力 的 影响 ,发 现 材料 粒度 越 小 ,切削 力 越 大 , 比 切削 
能 茵 高 。 这 是 由 于 随 着 材料 的 平均 粒 径 减 小 , 唱 界 
的 侣 响 增强 , 剪 切 力 和 天 切 力也 随 之 变 得 更 大 站 。 
Liu 等 中 对 不 同 晶 粒 尺寸 的 Al 5000 铝 合金 微 铣削 ， 
通过 细 化 材料 微观 结构 可 以 显著 改善 工 腔 的 表面 完 
整 性 。 残 余 应 力 是 评价 切削 完成 后 零件 表面 质量 的 
重要 指标 , 朱 黛 茹 "在 微 铣削 硬 铝 合金 LY12 实验 
基础 上 ,建立 了 基于 响应 曲面 法 和 回归 分 析 方 法 的 
微 铣削 残余 应 力 预测 模型 ,模型 预测 值 与 实验 测量 
值 吻合 ,该 模型 回归 显著 , 拟 合 度 达 到 97. 5% ,证 实 
模型 的 有 效 性 。 周 军 " 通 过 微 车 前 残余 应 力 实验 ， 
分 析 不 同 切削 深度 对 工件 圆周 方向 和 轴 向 方向 上 残 
余 应 力 的 影响 规律 ,发 现 工件 表面 两 个 方向 上 均 为 
残余 拉 应 力 ,并 且 拉 应 力 的 值 随 着 切削 深度 的 增 大 
而 增 大 。 路 疹 君 "以 镍 基 高 温 合 金 meconel718 为 研 
究 对 象 , 在 ABAQUS 里 建立 微 铣削 三 维 仿真 模型 ， 
研究 了 每 齿 进 给 量 对 进 给 方向 和 垂直 于 进 给 方向 残 


余 应 力 的 影响 规律 并 通过 实验 进行 验证 。 

材料 的 本 构 模 型 反映 的 是 应 力 与 温度 、 应 变 、 应 
变速 率 等 材料 状态 变量 的 映射 关系 。 在 传统 切削 研 
究 中 ,常用 的 本 构 模 型 无 法 解释 在 微 铣削 小 尺度 下 
的 尺寸 效应 和 最 小 切削 厚度 等 现象 。Liu 等 基 
于 泰勒 位 错 机 制 的 应 变 梯度 理论 ,修正 了 正 交 微 切 
削 铝 过 程 的 本 构 模 型 ,对 比 有 无 应 变 梯度 对 切削 力 
以 及 切削 温度 的 影响 ,发 现 修正 的 本 构 模 型 可 以 很 
好 地 表示 尺度 效应 。Wu 等 "将 应 变 梯 度 塑 性 理论 
应 用 到 45 钢 微 切 前 第 了 变形 区 模型 ,研究 主 剪 切 变 
形 区 的 有 效 流动 应 力 、 应 变 分 布 及 变化 规律 。 

残余 应 力 的 产生 是 一 个 非常 复杂 的 过 程 ,与 切 
削 加 工 过 程 中 的 热力 耦合 密切 相关 ,不 同 切削 参数 、 
刀具 加工 方 式 都 会 对 残余 应 力 产 生 较 大 的 影响 。 
上 述 学 者 在 对 残余 应 力 研 究 中 没有 将 工件 的 唱 粒 尺 
十 考虑 在 内 。 微 细 铣 削 中 ,刀具 刃 口 圆 绝 半径 、 切 削 
用 量 以 及 工件 材料 的 晶 粒 尺寸 都 处 于 微米 级 别 , 实 
际 上 的 切削 相当 于 在 不 同 唱 粒 中 进行 ,相当 于 安 观 
上 的 不 连续 切削 ,会 对 切削 过 程 中 塑性 变形 以 及 切 
削 热 产生 影响 ,进而 对 已 加 工 表面 残余 应 力 的 分 布 
和 大 小 产生 影响 。 

本 研究 将 通过 多 组 正 交 仿真 分 析 对 已 加 工 表 面 
残余 应 力 进 行 研究 ,建立 多 元 线性 回归 模型 并 进行 
显著 性 检验 ,通过 实验 对 预测 模型 准确 性 进行 验证 。 


1 仿真 模型 建立 及 验证 


1.1 微 铣削 过 程 模型 简化 


微 铣削 加 工 是 一 个 断 续 切削 过 程 , 加 工 过 程 中 切 
削 深 度 时 刻 随 着 刀 夷 的 运动 而 变化 。 实 际 三 维 铣削 
仿真 由 于 切削 过 程 ,接触 条 件 的 设置 复杂 ,其 仿真 速度 
慢 ,计算 结果 不 容易 收敛 ,因此 在 保证 仿真 结果 准确 的 
前 提 下 ,可 以 通过 合理 的 假设 对 仿真 模型 进行 简化 。 

基于 平面 应 变 状 态 假 设 条 件 :中 切削 宽度 大 于 
等 于 5 倍 切 前 深度 ;@ 塑 性 变形 仅 发 生 在 与 切削 丸 
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垂直 的 平面 。 当 满足 以 上 两 个 假设 条 件 时 ,可 以 根 
据 弹 塑性 力学 理论 将 铣削 加 工 简化 为 平面 应 变 问 
题 ,从 而 可 以 利用 二 维 正 交 切削 模型 进行 仿真 分 析 。 

由 于 切削 深 度 非 常 小 ,可 以 忽略 铣 刀 切削 刃 的 
螺旋 角 , 因 此 可 以 将 图 1(a) 所 示 的 实际 微 铣削 过 程 
简化 为 图 1(b) 所 示 的 平面 模型 ,即将 三 维 微 铣削 过 
程 中 的 恒定 切 前 深度 转变 为 平面 模型 中 的 工件 厚度 。 
加 工 变形 区 的 太 才 (微米 级 和 亚 微米 级 ) 远 小 于 微 铣 刀 
直径 ( 几 百 微米 ) ,根据 图 1(e) 所 示 的 铣削 过 程 工 件 变 
形 区 局 部 放大 图 ,其 类 似 于 图 1(d) 的 正 交 切削 过 程 ， 
此 时 切削 深度 相当 于 微 铣削 过 程 中 每 齿 进 给 量 ” 。 


二 
进 给 方向 
(a) 微 铣削 实际 过 程 


(b) 微 铣削 简化 平面 模型 


a (c) 局 部 放大 医 


(d) 正 交 切削 模型 

.二 图 1 微 铣削 过 程 简 化 

mm Fig.1 Micro milling process simplification 

@ 图 1 中 :r 为 刀具 刃 口 国 弧 半径 jj 为 每 齿 进 给 
量 ; 刀 ,为 正 应 力 ; 挛 为 摩擦 力 。 


1.2 多 晶体 几何 模型 建立 


基于 Voronoi 图 的 结构 原理 "1 ,通过 MATLAB 
生成 随机 的 唱 粒 结构 ,并 将 其 导入 有 限 元 软件 
ABAQUS 生成 工件 几何 模型 。 图 2 为 包含 18 个 唱 
粒 的 多 晶体 几何 模型 ,其 平均 唱 粒 尺寸 为 33.3 nm， 
其 完整 的 程序 流程 图 如 图 3 所 示 。 
1.3 微 铣削 物理 模型 
1.3.1 本 构 模 型 

材料 的 本 构 模 型 描述 的 是 应 力 与 温度 .应变 应 
变速 率 等 材料 状态 变量 的 变化 情况 ,可 以 反映 变形 
过 程 中 材料 本 质 的 变化 ,对 有 限 元 仿真 效果 及 成 败 
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起 着 决定 性 的 作用 。 
将 泰勒 位 错 理 论 和 应 变 梯度 理论 引入 ,在 JC 本 
构 模 型 的 基础 上 进行 修正 ,推导 如 下 。 
泰勒 位 错 理论 描述 剪 切 流动 应 力 7 与 位 错 密度 
Pp 的 关系 为 
7 =abG yp =abG Vp, (1) 
0 = /7 (2) 
统计 存储 位 错 p. 与 材料 单 轴 拉 伸 参 考 流动 应 力 
0 的 关系 为 
au = kabG ps = 0 (3) 
微 切 前 中 由 于 尺寸 效应 的 影响 ,材料 的 剪 切 流 
动 应 力 由 统计 存储 位 错 p. 和 几何 必需 位 错 pe 共 同 决 
定 , 于 是 可 以 将 式 (1) 变 换 为 式 (4) , 即 
7 =abG Vp =abG Vp. +pe (4) 


图 2 多 晶体 几何 模型 


Fig.2 Polycrystalline geometric model 


生成 晶 粒 数 的 随机 点 


利用 MPT 工 具 箱 ， 生 成 
Voronoi 图 的 基本 信息 


Matlab 
程序 
部 分 


到 的 数据 按 点 、 线 、 
顺序 输出 到 txt 文 本 


读 取 txt 文 本 中 数据 ， 存 放 
于 多 级 列表 


Python 三 调用 列表 妆 据 ， 利用 WirePloyine 
神 序 站 令 将 点 生成 线 ， 利 用 AddFace 
部 分 命令 将 线 生 成 面 


图 3 多 晶体 几何 模型 建立 的 流程 图 
Fig.3 Flow chart of the establishment of 


a polycrystalline geometric model 


几何 必需 位 错 可 由 应 厚 探 模型 


应 变 梯度 7 表示 , 即 
Be (5) 
联 立 可 以 得 到 考虑 应 变 梯度 的 流动 应 力 , 即 
og=0k 1+2 (0) bm (6) 


Li 等 “研究 得 出 当 晶体 发 生 塑 性 变形 时 ,其 应 变 
2 统计 存储 位 错 P. 和 品 粒 太 十 4 之 间 的 近似 关系 为 


Fog 
ps = 714 (7) 


联 立 以 上 各 式 可 以 得 到 适用 于 多 品 体 材料 切削 
的 本 构 方程 , 即 


aicd\ 
2 1+2 | Oe ) ji (8) 


式 中 :uv 为 JC 模型 下 的 等 效 流动 应 力 ;a 为 实验 党 
数 * s 昌 Qa =0.5;b 为 Burgers 矢量 ;6 为 剪 切 模 量 ,C = 
86. 003 GPa;”m 为 应 变 梯度 , 且 m =0. 38'™。 

回 式 (8) 为 在 JC 模型 的 基础 上 ,基于 细 观 力学 应 
区 稀 度 理 ; 企 和 位 错 理论 建立 的 适用 于 多 晶体 材料 微 
十 前 过 程 中 材料 本 构 模 型 ,其 反映 了 应 力 与 应 变 、 应 
误 素 温度 的 关系 ,并 且 描 述 了 微 切 前 过 程 中 尺寸 效 
应 等 特有 现象。 根据 VUMAT 子 程序 二 次 开发 计算 
和， 实施 的 总 体 流程 图 如 图 4 所 示 。 


人 
调用 VUMAT 


输入 材料 参数 (JC 
参数 及 应 变 梯度 参 
数 为 s 


前 应 
变 增 量 


计算 弹性 试探 应 力 


计算 静水 应 力 、 偏 应 
力 及 Mises 等 效应 力 


妈 4” 微 铣 前 本 构 方 程 子 程序 开发 流程 区 


Fig.4 Flow chart of development of micro-milling 


constitutive equation subroutine 


1.3.2 摩擦 模型 
金属 切削 过 程 中 ,在 第 开 变 形 区 和 第 亚 变 形 区 
都 存在 摩擦 接触 的 关系 ,摩擦 的 作用 会 导致 刀具 强 
度 和 寿命 降低 、 对 切削 过 程 中 切削 力 和 切削 温度 产 
生 较 大 的 影响 ,从 而 对 工件 的 表面 加 工 质 量 和 加 工 
精度 产生 严重 的 影响 。 

在 金属 切削 有 限 元 仿真 软件 ABAQUS 里 通常 
使 用 的 摩擦 模型 为 修正 的 库伦 摩擦 模型 ,该 模型 可 
以 很 好 的 体现 黏 结 区 和 滑 移 区 刀具 和 工件 的 摩擦 接 
触 情况 ,其 数学 表达 式 为 


“(滑动 区 ) 


LO Tan <Th 
fF 9 
ee ,Tn 宇 Tnax( 医 结 区 ) 
式 中 : 7 摩擦 应 力 ;o, 为 正 应 力 ;7, 为 工件 材料 的 


风 服 应 方 ;4 为 摩擦 系数 。 
1.3.3 切 悄 分 离 准 则 
金属 切削 仿真 模拟 分 析 中 切 悄 和 工件 的 分 离 准 
则 主要 包括 几何 分 离 和 物理 分 离 。 物 理 分 离 准 则 中 
的 JC 断裂 准则 可 以 同时 反映 出 被 加 工 材料 的 物理 
性 质 和 力学 性 质 ,使 有 限 元 仿真 模拟 更 接近 于 真实 
情况 ,其 表达 式 为 
六 全 [| (10) 
Er 

式 中 :A 5 为 等 效 塑 性 应 变 增 量 ;2 为 等 效 破坏 应 变 。 

si 的 表达 式 如 式 (11) 所 示 , 失 效 参 数 如 表 2 
所 示 。 


E 了 一 了 
过 [a + dexp( d; | +dan( 2 )]l! | | 


(11) 
式 中 :p 为 静水 压 应 力 ;g 为 偏 应 力 ;di ~ ds 为 失效 参 
数 ,如 表 1 所 示 。 
表 1 GH4169 合金 的 失效 参数 
Tab.1 Failure parameters of GH4169 alloy''® 
失效 参数 di d, ds da ds 
数值 0.239 0.456 -0.3 0.07 2.5 
1.3.4 热传导 模型 


由 于 塑性 变形 和 刀 - 眉 摩擦 会 产生 大 量 切 前 热 ， 
但 是 加 工时 间 非 常 短 ,通常 只 考虑 工件 和 刀具 之 间 
的 热传导 ,而 忽略 工件 与 周围 介质 的 热 对 流 及 工件 
与 刀具 表面 的 热 辐射 影响 。 二 维 切 前 仿真 中 热传导 
方程 为 


oT 97 a7 97 
+a] + =p.C,l| 1 (12) 
式 中 :A 为 热 导 率 ; C0, 为 比热容 ;uv 为 动态 热源 在 
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和 了 方向 上 分 支 ;9 为 单位 体积 内 总 热 生产 率 ;p。 为 
工件 的 密度 。 
1.3.5 边界 条 件 
在 进行 微 铣削 残余 应 力 仿真 时 ,将 工件 底部 和 
左 侧面 完全 国定 ,给 刀具 参考 点 施加 沿 x 负 方 向 初 
台 速 度 ,其 大 小 为 125. 6 mmin。 同 时 为 了 更 加 契合 
实际 铣削 过 程 ,需要 对 工件 和 刀具 设置 温度 边界 条 
件 ,将 其 所 有 节点 的 初始 温度 定义 为 室温 , 取 25 %C。 


1.4 仿真 模型 验证 


为 了 验证 所 建立 的 不 同 晶 粒度 的 多 晶体 仿真 模 
型 是 否 准确 ,将 基于 Voronoi 图 建立 的 多 晶体 几何 模 
型 进行 切削 仿真 ,将 其 仿真 得 到 的 切削 力 与 微 铣削 
实验 中 测 得 真实 的 切 前 力 进行 比较 ,得 到 其 相对 误 
差 3s 本 次 实验 采用 汉 川 机 床 公司 生产 的 三 坐标 轴 立 
控 铣 床 ,型 号 为 XK714D。 加 工 的 材料 为 经 过 
热 您 理 后 平均 晶 粒 尺寸 为 20 pm 的 镍 基 高 温 合金 
CEB169 板材 ,利用 Kistler 测 力 仪 采集 4 组 微 铣 档 切 
前 放 。 微 铣 槽 实验 参数 如 表 2 所 示 。 刀 具 采 用 的 株 
洲 便 质 合金 厂 生产 的 4 刃 平底 硬 质 合金 立 铣 刀 ,型 
E 0s ,法 前 角 y, = 10。, 后 角 w, =8°, 晶 
旋 莉 5, = 45°, 涂 层 材 料 为 立方 氮 化 硼 ,不 采用 冷 
却 液 。 
表 2 仿真 模型 验证 微 铣 槽 实验 参数 


ac Tab.2 Simulation model to verify the experimental 


parameters of micro-milling slots 


和 。 和 和 于 。 每 寺 进 给 量 / 切削 深度 / 
四 Pp™ HRm pm 
O 4 000 9 时 
2 5 000 7 
3 6 000 5 人 
4 7 000 3 加 


仿真 过 程 中 切削 深度 是 保持 在 名 义 上 每 齿 进 给 
量 不 变 的 ,仿真 和 实验 得 到 的 式 和 了 切削 力 代 表 含 
义 不 同 ,不 能 直接 形成 对 比 ,因此 本 研究 对 X 和 了 方 
向 的 合力 进行 验证 。 其 切削 力 对 比 及 误差 如 表 3 
所 示 。 

由 表 3 可 知 , 基 于 修正 的 JC 本 构 模 型 和 常规 JC 
本 构 模 型 的 仿真 与 实验 切削 力 相对 误差 最 小 的 是 第 1 
组 ,其 值 分 别 为 6.32% 和 26.07% ,误差 最 大 是 第 4 组 
为 16.18% 和 38.58% ,平均 误差 分 别 为 10.92% 和 
31. 88% 。 基 于 修正 的 JC 本 构 模 型 的 切削 力 误差 明 
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显 小 于 常规 JC 本 构 模 型 , 且 随 着 每 从 进 给 量 的 减 
小 ,基于 常规 的 J 本 构 模 型 的 切 前 力 误差 呈现 增 大 
趋势 ,这 主要 是 由 于 微 切 前 中 尺寸 效应 导致 的 。 而 
基于 修正 的 JC 本 构 模 型 的 切 前 力 误 差 一 直 很 小 ,可 
以 很 好 地 解释 尺寸 效应 现象 ,从 而 验证 本 研究 所 建 
立 仿真 模型 的 准确 性 。 

表 3 仿真 模型 验证 误差 


Tab.3 Simulation model verification errors 


修正 模型 JC 模型 
,。 ”实验 值 / 误差 / 误差/ 

编号 仿真 值 / 仿真 / 

N % % 
N N 

1 21.98 23.37 6.32 16.25 26.07 
27.82 31.38 12.79 36.10 29.76 
3 17.44 18.91 8.41 235:22 33. 14 
4 23.05 26.78 16.18 31.94 38.58 


2 残余 应 力 仿真 过 程 


残余 应 力 的 有 限 元 仿真 过 程 一 般 分 为 4 个 阶 
段 :切削 加 工 阶段 、 退 刀 阶 段 、 约 东 转 换 以 及 冷却 
阶段 。 
S,Mises 


(Avg:75%) 
二 1.093e+03 


图 5 切削 加 工 阶段 完成 后 Mises 应 力 云图 


Fig.3 Mises stress cloud after completion of the cutting process 


S,Mises 
(Avg:75%) 
+8.197e+02 


图 6 退 刀 和 转化 约束 后 Mises 应 力 云图 
Fig.6 Mises stress cloud after tool withdrawal and 
transformation constraints 
为 了 获取 已 加 工 表面 残余 应 力 值 ,在 水 平方 向 
上 等 间隔 取 10 个 节点 ,然后 在 ABAQUS 后 人 处理 中 提 


取 其 对 应 点 的 数值 ,将 所 提取 出 数值 取 平 均值 得 到 
已 加 工 表面 的 残余 应 力 值 。 
S,Mises 


(Avg:75%) 
+7.625e+02 


图 7 自然 冷却 后 Mises 应 力 云 图 


Fig.7 Mises stress cloud after natural cooling 


S,S11 

(Avg:75%) 
+6.000e+02 
+$5.077e+02 


图 8 最 终 残余 应 力 云图 


Fig.8 Final residual stress cloud 


| UR Ue te a 
‘oO 
OI 一 上 中 一 上 一 
上 OOOO%WPEUD 
Be 小 上 上书 一 局 上 个 
Ooaqcaqaaaqaqaacm 
士 十 十 十 十 十 十 十 于 守 字 
OOCOCOOOOOOOO 
DO 


3 已 残 余 应 力 预测 模型 


2 


~ 
3 残余 应 力 正 交 仿真 


四 本 研究 主要 研究 刀具 刃 口 圆 弧 半径 ,切削 深度 

粒 尺 寸 对 微 铣削 高 温 合金 GH4169 已 加 工 表面 

应 力 的 影响 规律 。 设 计 三 因素 四 水 平 正 交 仿 真 
实验 ， 其 因素 水 平 表 如 表 4 所 示 。 根 据 因素 水 平 
表 ， 选 用 标准 正 交 实验 表 L16 (45) 设计 正 交 实验 
仿真 方案 来 研究 各 因素 对 微细 铣削 高 温 合金 
CH4169 表面 残余 应 力 的 影响 ， 正 交 表 及 其 仿真 结 
果 如 表 5 所 示 。 

表 4 正 交 仿真 因素 水 平 表 


Tab.4 Level table of orthogonal simulation factors 


A 晶 粒 尺 寸 / B 刃 口 半径 / C 切削 深度 / 
水 平 
Hm Hm Hm 
1 10 5 4 
2 20 10 9 
3 33.3 20 15 
4 40 30 25 
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表 5 正 交 仿真 方案 及 结果 


Tab.5 Orthogonal simulation scheme and results 


编号 A B C 水 平 组 合 OZMPa 
1 1 1 1 AlBI1Cl1 258.74 
1 2 2 AlB2C2 214.32 
3 1 3 3 AlB3C3 275.83 
4 1 4 4 AlB4C4 231.98 
| 2 1 2 A2B1C2 190. 18 
6 2 2 1 A2B2C1 248.59 
7 2 3 4 A2B3C4 200.37 
8 2 4 3 A2B4C3 180.63 
9 3 1 3 A3B1C3 189.82 
10 3 2 4 A3B2C4 170.13 
11 3 3 1 A3B3C1 198.25 
12 3 4 2 A3B4C2 190.49 
13 4 1 4 A4B1C4 170.71 
14 4 2 3 A4B2C3 194.33 
15 4 3 2 A4B3C2 195.45 
16 4 4 1 A4B4C1 223. 69 


3.2 ”残余 应 力 预测 模型 建立 


假设 已 加 工 表面 残余 应 力 与 晶 粒 尺寸 .切削 深 
度 以 及 刀具 为 口 圆 弧 半径 之 间 存 在 复杂 的 指数 关 
系 ,通过 统计 学 的 方法 建立 残余 应 力 的 正 交 回归 的 
预测 模型 , 即 


ai = Cd ar (13) 
式 中 :0 ,为 残余 应 力 ;C 为 常数 ,数值 大 小 由 加 工 材 
料 和 切削 条 件 决定 ;81 ,5, ,5; 为 待定 系数 。 
通过 变化 得 到 其 相应 的 线性 回归 方程 为 
y=bo +bx) 十 DOX， + Osx (14) 
在 回归 方程 中 加 上 一 项 实验 误差 , 则 16 组 仿真 
结果 建立 的 多 元 线性 回归 方程 可 以 用 和 矩阵 形式 表示 为 


Y=XB +e， 
Yi 1 Xl X12 X13 21 
J 1 MX21 N22 N23 22 
Y= ,在 三 ,2 二 和 

yi16 1 X161 Ni62 Ni63 £16 

Bo 

bi 

Bb; 

pb; 


参数 B 可 通过 最 小 二 乘法 进行 估计 , 则 回归 方 
程 的 估计 量 为 
y=bo +bix, + bx + baxs (15) 
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式 中 bo ,b1,b, ,bs 分 别 对 应 参数 Bo ,GD ,B; 的 最 小 
二 乘 估计 , 则 
bp=(X'X)-1XY (16) 
将 向 量 半 和 了 代入 式 (18) ,可 得 
2.1726 
-0.055 7 
0. 105 5 
0. 112 3 
则 根据 正 交 仿真 结果 所 得 到 的 数据 ,得 到 残余 
应 力 的 预测 模型 , 即 
Cr =148. 798 9d "a WE (17) 


3.3 ”预测 模型 显著 性 检验 


。_ 在 建立 残余 应 力 预 测 模型 之 前 ,随机 变量 y 与 
oo va swi 之 间 的 关系 不 确定 ,只 能 通过 假设 的 方 
法 来 进行 判定 。 因 此 当 求 出 预测 模型 之 后 ,需要 对 
模 型 的 显著 性 进行 检验 以 此 来 判断 预测 模型 拟 合 的 
妤 环 程 度 。 

已 为 了 进行 统计 检验 ,将 回归 模型 的 总 偏差 平方 
人 表示 为 


CO Sr =SA+TSF (18) 
式 计 :S\ 为 回归 偏差 平方 和 ;Se 为 剩余 平方 和 ,可 以 
表 丰 为 


© Sr = 2, Cy; -7 ， 
三 Sa 一 > (7 7 


中 采用 检验。 假设 包 :B, =0,B, =0,B; =0, 可 
采用 式 (19) 所 示 的 统计 量 , 即 
SA/s 
Se/(n-s-1) 
本 研究 实验 组 数 n=16, 变 量 个 数 s=3, 分 别 计 
算式 (18) 和 式 (19), 显 著 因 子 取 0.01, 将 计算 结果 
列 于 表 6。 查 分 布 表 得 Fo0(3,12) =5.95 < 下 = 
25.43 ,拒绝 原 假设 ,因此 本 研究 所 建立 的 残余 应 力 
预测 模型 是 高 度 显 著 的 。 
表 6 残余 应 力 方差 分 析 表 


Tab.6 Variance analysis table of residual stress 


F ~F(Si,n-s-1) (19) 


方差 来 源 平方 和 由 度 方差 了 值 


SA 17.54 3 5.85 25.43 
SE 2.76 12 0.23 
Si 20.30 15 四 
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4 预测 模型 实验 验证 


当 合金 成 分 一 定时 ,热处理 工艺 就 成 为 改变 合 
金 组 织 的 重要 手段 。 不 同 的 热处理 方案 ,可 以 控制 
合金 的 晶 粒 尺寸 ,析出 强化 相 的 数量 、 形 貌 , 进 而 可 
以 根据 实际 需求 , 使 合金 具有 不 同 的 力学 性 能 。 
GH4169 合金 最 常用 的 热处理 工艺 需要 经 历 固 浴 处 
理 和 双 级 时 效 处 理 两 个 阶段 。 其 中 国 溶 处 理 有 两 个 
目的 ,其 一 是 溶解 材料 中 基本 组 成 相 和 杂质 物 , 从 而 
形成 一 种 均匀 的 过 饱和 固 洲 体 ,其 二 是 为 了 改变 唱 
粒 尺 寸 。 本 研究 确定 其 热处理 方案 如 表 7 所 示 , 其 
最 终 测 得 平均 晶 粒 尺寸 分 别 为 10 pm 和 40 pm。 

表 7 热处理 工艺 


Tab.7 Heat treatment process 


固 溶 处理 


1 985%C xlh, 空 准 


时 效 处 理 
720 % xlh 炉 冷 +620 CC x8h 空冷 


2 1035% xlh, 空 冷 ” 720% xlh 炉 冷 +620% x8h 空冷 


4.1 实验 装备 


如 图 9 所 示 , 实 验 所 用 机 床 为 日 本 NanoWare 公 
司 生产 的 精密 数控 微 铣床 , 型 号 为 MISSR ,工件 的 
安装 及 预 加 工 如 图 10(a) ~ 10(b) 所 示 。 检 测 设备 
是 日 本 Rigaku International Corpratio 公司 生产 的 X 
射线 衍射 仪 ,其 型 号 为 DVMAX2500HB +ZPC, 其 中 
XRD 实验 的 参数 设置 : 靶 材 为 Cu、 准 直 管 直径 为 
0.5、 镜 面 为 311、 衡 射 角 为 20° ~ 90°、 弹 性 模 量 为 
206 GPa 和 泊 松 比 为 0.3。 


图 9 MTS5R 精密 数控 微 铣床 
Fig.9 MTSSR precision CNC micro-milling machine 
实验 工件 材料 为 钊 基 高 温 合金 CH4169 ,工件 形 


状 为 B21 mm x7 mm 圆 饼 状 ,中 间 为 G5 mm 的 通 
孔 , 如 图 12(e) 所 示 。 刀 有 具 采用 定制 的 日 本 NS Tool 
公司 生产 的 4 为 平底 硬 质 合金 立 钳 刀 , 型 号 为 
MXH225 , 刃 口 圆 弧 半径 分 别 为 10 pm 和 20 km, 其 


他 主要 参数 : 涂 层 为 TiAIN 、 刃 径 为 1 mm、 刃 长 为 
5 mm \ 刀 长 为 45 mm 、 锥 柄 角 为 9 "和 螺旋 角 为 25"。 


们 7 “i 时 
(4) 


(b) 预 加 工 


图 10 工件 安装 及 预 加 工 


Fig. 10 Workpiece installation and pre-processing 


(a) 安装 


4.2 实验 方案 


实验 切削 速度 与 仿真 中 一 样 , 均 为 125.6 m/min， 

其 唱 粒 尺寸 、 刃 口 半 径 和 每 齿 进 给 量 (由 1.1 节 切 削 

仿 蛋 模型 简化 中 切削 深度 相当 于 微 铣削 过 程 中 单 齿 

给 量 ) 参 数 的 选择 在 正 交 仿 真 选 取 的 各 因素 的 范围 
光 随 机 选取 4 组 参数 进行 实验 验证 ,如 表 8 所 示 。 
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表 8 残余 应 力 模型 验证 实验 表 


Tab.8 Residual stress model verification experiment table 
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水 平 ” 晶 粒 尺寸 /jm ” 刃 口 半径 /pm 每 从 进 给 量 /hm 
1 10 10 5 
2 10 10 10 
3 40 20 15 
4 40 20 20 
4.3 实验 结果 
将 样品 进行 全 谱 扫 描 , 确 定 合适 的 衍射 峰 , 然后 


按照 提前 设置 好 的 4 个 侦 转 角度 慢 速 扫描 选择 衍射 
峰 , 测 出 图 谱 后 采用 JADE 软件 中 的 定 峰 方法 得 到 
20 值 ,在 ORIGIN 里 画 出 并 拟 合 得 到 关系 曲线 ,得 到 
曲线 的 斜率 ,进而 可 以 求 出 残余 应 力 值 。 

表 8 所 示 切 削 参数 加 工 后 的 表面 经 过 XRD 扫 
描 后 图 谱 如 图 11 所 示 , 其 测量 结果 中 , 出 现 衍射 峰 
的 角度 的 分 布 情况 如 表 9 所 示 。 


©O 7 500 
国 
号 10 000 6 000 
全 8 000 4 500 
萎 
ON 6 000 绰 
加 图 ad 
~> 村 3 000 
a 4 000 
>< 1 500 
© 2 000 
9 
re 20 30 40 5 60 70 80 90 20 30 40 50 60 70 80 90 
20 20 
© (a) 第 1 组 (b) 第 2 组 
9 000 
8 000 
7 500 
6 000 
6 000 
Em 
量 4500 = 4 000 
3 000 
200 
1 500 
0 0 
20 30 40 5 60 70 80 90 20 30 40 50 60 70 80 90 100 
20 20 
(c) 第 3 组 (d) 第 4 组 
图 11 实验 对 照 组 XRD 衍射 图 谱 


Fig. 11 


XRD diffraction pattern of the experimental control group 
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表 9 产生 衍射 峰 的 20 值 
Tab.9 20 value of diffraction peak 


实验 号 产生 衍射 峰 的 位 置 26 
1 47.28 54.60 75.34 86.42 
2 43.52 50.68 73.90 87.96 
3 43.66 50.54 74.44 87.82 
4 46.24 56.59 76.26 87.58 


siny siny 


siny 
(d) 第 4 组 
图 12 ”实验 对 照 组 20-sinyp 拟 合 曲线 
Fig.12 ”Experimental control group fitting curve 
表 10 残余 应 力 模型 误差 


Tab. 10 Residual stress model errors 


预测 值 /MPa 实验 值 /MPa 误差 /% 
200.835 7 210.0942 4.61 
Ee 216. 121 7 231.9707 7.35 
243.8159 238.5007 2.18 
251.328 1 262.048 1 4.29 


表 10 为 通过 正 交 仿 真 建立 的 回归 模型 与 XRD 
实验 测 得 的 残余 应 力 值 的 误差 , 其 误差 分 布 在 
2.18% ~7.35% ,不 超过 10% ,因此 可 以 说 明 残 余 应 
力 的 预测 模型 准确 性 较 高 ,对 于 后 期 实际 加 工 过 程 
中 残余 应 力 的 预测 有 一 定 的 指导 作用 。 


5 结 论 


1 ) 根 据 平面 应 变 状态 假设 条 件 , 将 其 简化 为 二 
维 正 交 切 前 模型 ,基于 Voronoi 图 的 结构 原理 建立 多 
晶体 GH4169 二 维 几 何 模 型 。 将 同 条 件 JC 模型 \ 修 
正 后 JC 模型 得 到 的 仿真 切削 力 与 实验 测 得 值 对 比 
并 计算 相对 误差 ,误差 最 小 的 分 别 为 26. 07% 和 
6. 32% ,误差 最 大 是 38. 58% 和 16. 18% ,平均 误差 
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分 别 为 31.88% 和 10.92% 。 基 于 修正 的 JC 模型 的 
切削 力 误差 小 于 JC 模型 且 一 直 很 小 ,可 以 很 好 地 解 
释 尺 寸 效 应 ,从 而 验证 了 本 研究 所 建立 的 仿真 模型 
的 准确 性 。 

2) 运 用 多 元 线性 回归 理论 对 残余 应 力 与 晶 粒 尺 
寸 . 切 削 深 度 、 刀 具 为 口 圆 弧 半 径 的 数学 关系 进行 了 
预测 ,并 对 模型 的 显著 性 进行 检验 ,结果 表明 模型 是 
高 度 显 著 的 。 

3 ) 对 高 温 合金 GH4169 进行 不 同 热处理 工艺 ， 
得 到 不 同 晶 粒 尺 寸 , 设 置 与 仿真 参数 一 致 的 微 铣削 
实验 ,运用 XRD 测量 表面 残余 应 力 ,将 实验 测量 的 
残余 应 力 与 模型 预测 的 进行 误差 分 析 , 其 误差 分 布 
在 2.18% ~7.35% ,不 超过 10% ,因此 可 以 说 明 残 
余 应 力 的 预测 模型 准确 性 较 高 ,对 于 后 期 实际 加 工 
过 程 中 残余 应 力 的 预测 有 一 定 的 指导 作用 。 
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